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Abstract. We experimentally characterize the non-equilibrium, room-temperature 
magnetization dynamics of a spin chain evolving under an effective double-quantum 
Hamiltonian. We show that the Liouville space operators corresponding to the 
magnetization and the two-spin correlations evolve 90 degrees out of phase with each 
other, and drive the transport dynamics. For a ncarcst-neighbor-coupled A-spin 
chain, the dynamics are found to be restricted to a Liouville operator space whose 
dimension scales only as A 2 , leading to a slow growth of multi-spin correlations. 
Even though long-range couplings are present in the real system, we find excellent 
agreement between the analytical predictions and our experimental results, confirming 
that leakage out of the restricted Liouville space is slow on the timescales investigated. 
Our results indicate that the group velocity of the magnetization is 6.04 ± 0.38 /xm/s, 
corresponding to a coherent transport over N m 26 spins on the experimental timescale. 
As the double-quantum Hamiltonian is related to the standard one-dimensional XX 
Hamiltonian by a similarity transform, our results can be directly extended to XX 
quantum spin chains, which have been extensively studied in the context of both 
quantum magnetism and quantum information processing. 



PACS numbers: 03.67.Hk, 75.10.Pq, 76.90.+d 
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1. Introduction 

Solid-state spin systems provide an attractive test-bed to study both equilibrium 
and non-equilibrium quantum many-body dynamics, and have recently emerged as 
a promising platform for quantum simulation [1, 2]. One-dimensional (ID) spin 
systems, in particular, are of special interest as they connect to an important class 
of problems in condensed-matter physics [3], and have been suggested as quantum 
wires to coherently transfer quantum information across distant nodes in a quantum 
computer and distributed quantum architectures [4, 5]. Transport properties, and 
magnetization transport in particular, have been extensively investigated theoretically 
by both the condensed-matter community and more recently in the context of spintronics 
and emerging nano-device applications [6]. 

The class of ID spin-1/2 XY Hamiltonians [7, 8, 9], which are exactly solvable via 
a Jordan- Wigner mapping onto a system of non-interacting spinless fermions, play an 
archetypal role in condensed-matter physics. In this case, a local magnetic disturbance 
is known to propagate down the chain scatter-free with a constant velocity {ballistic 
transport), rather than diffusively spreading from the site of the disturbance and 
eventually decaying {diffusive transport). From a quantum communication perspective, 
the mapping to free fermions has proved crucial to also allow a quantum state to be 
transported down the chain — which thus acts as a quantum conduit or channel (see 
for example [10] for a recent review). While at zero temperature, within linear response 
theory, integrable quantum models are typically associated with ballistic transport, a 
coexistence with and/or crossover to diffusive behavior may be possible more generally, 
for instance in the presence of non-local conserved quantities [11] or of couplings to an 
environment [12, 13, 14]. Despite significant progress, a satisfactory understanding of 
the conditions leading to ballistic versus diffusive transport is as yet lacking, with a 
number of unresolved questions remaining, in particular, in relation to the impact of 
finite and infinite temperatures [15], the role of non-integrability, and its interplay with 
frustration [16, 17]. As a result, experimentally characterizing the transport properties 
and physical mechanisms in low-dimensional spin systems remains important from both 
fundamental and applied standpoints. 

The long coherence times afforded by nuclear spins and the ability to access a large 
Hilbert space, in conjunction with the superb level of control available over spin degrees 
of freedom, make solid-state nuclear magnetic resonance (NMR) an excellent setting 
for exploring the coherent dynamics of a (nearly) isolated quantum many-body system 
[1, 18, 19, 20] as well as the statistical physics of equilibrating spin systems [21, 22, 23]. In 
this work, we employ solid-state NMR techniques to characterize the room-temperature 
magnetization dynamics of a (quasi) ID spin system during coherent evolution under an 
effective Double Quantum (DQ) Hamiltonian, which is directly related to the isotropic 
XY (XX henceforth) Hamiltonian by a similarity transformation. In particular, we 
find that the experimental results are in very good agreement with those predicted by 
the free-fermionic solutions, indicating that integrability-breaking perturbations (due to 
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longer range couplings) have negligible effect on the timescale of the experiments, and 
enabling us to calculate the transport velocity of the magnetization. This is in marked 
contrast to the diffusive behavior observed under the dipolar Hamiltonian in previous 
experiments on 3D spin networks [24, 25], where spin diffusion arises from unitary 
dynamics under the high-field (secular) dipolar Hamiltonian. While experimentally 
the high degree of isolation from the surrounding environment during evolution was 
demonstrated by the observation of "polarization echoes" upon reversing the sign of the 
dipolar Hamiltonian [26], quantum chaoticity was explicitly invoked theoretically as a 
mechanism for diffusion [27, 28]. 

Beside providing additional insight into the mechanisms underlying coherent 
transport in isolated ID many-body quantum systems, in the context of quantum 
information transport it is important to stress that the nuclear spin chains we study here 
are initially in a highly mixed quantum state. The importance of relaxing initialization 
constraints is being increasingly appreciated within the quantum communication 
community (see in particular [ ] and references therein). We thus expect this study 
to also be of direct relevance to a number of other quantum platforms where mixed- 
state spin chains are naturally encountered, such as phosphorus defects in silicon 
nanowires [ 0], quantum dots [31, 32], molecular semiconductors [ ] and solid-state 
defects in diamond or silicon carbide [34, 35]. 

2. Experimental Methods and Results 

The 19 F spins in a crystal of fluorapatite (FAp - Ca5(P0 4 ) 3 F) have long been used 
to experimentally approximate a nearest-neighbor coupled ID spin system (see for 
example [36, 37, 38, 39, 40]). In a 3D lattice of dipolar-coupled nuclear spins, every 
pair of spins is coupled with an interaction strength (between spin j and spin €) 
dj£ = (/i /167r)(7 2 ^/r| £ )(l — 3 cos 2 Ojg) [41], where 7 is the gyromagnetic ratio of fluorine, 
rj£ is the distance between nucleus j and £, and 9jt is the angle between Tji and the 
z-axis (along which the external magnetic field is applied). The geometry of the spin 
system is reflected in the distribution of the djg couplings. In FAp, the 19 F nuclei form 
linear chains along the c-axis, each one surrounded by six other chains. The distance 
between two intra-chain 19 F nuclei is r = 3.442 A whereas the distance between two 
cross-chain 19 F nuclei is R = 9.367 A. The largest ratio between the strongest intra- and 
cross- chain couplings (~ 40) is obtained when the crystalline c-axis is oriented parallel 
to the external field. In this orientation, the fluorine spins may be treated as a collection 
of many identical ID chains with only nearest-neighbor (NN) couplings, with a coupling 
strength d = (yU /8vr)7 2 /r 3 = 8.17 x 10 3 rad/s. 

In this NN approximation, the high-field secular dipolar Hamiltonian of a single 
chain is given by "Hoip = \ ^Yl!i=\ & {^ a i a i+i ~ ®i ' where o~ a (a = x,y, z) are the 

Pauli operators. Starting from %Di P and using suitable multiple-pulse sequences [43, 44], 
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we can experimentally implement an effective DQ Hamiltonian, given by 

N-l 

i 

Formally, the DQ Hamiltonian is related to the standard XX Hamiltonian by the 
similarity transformation U^q = exp(— in/2 ^\ <r*), where the sum is restricted to 
either even or odd spins. While I z oc ^ P z is conserved under the XX Hamiltonian, the 
corresponding conserved quantity under the DQ Hamiltonian is I z oc ^dq^^dq^- This 
results in the inversion of sign of the local magnetic disturbance at every alternate site 
as it moves down the chain, as will be seen later (see Figure 6). 

In order to implement the DQ Hamiltonian, we used a standard 8-pulse sequence 
applied on-resonance to the 19 F Larmor frequency [43]. This 8-pulse sequence, S = 
C-C-C-C, may be understood in terms of a simpler 2-pulse cycle (C, and its time-reversed 
version C), which also simulates the DQ Hamiltonian. The primitive pulse cycle is given 
by C = [|XA'Xy], where A' = 2A + w, A is the delay between pulses and w is the 
width the n/2 pulse (denoted by X), applied about the x-axis. The dynamics in the 
presence of the pulse sequence can be expressed in terms of a time-independent effective 
Hamiltonian Hoq, 

U x MQ (t) = T exp f-i J\u dip + H* f {s)\d^ = e-*W (2) 

where T denotes time-ordering operator, h — 1, and H^At) is the time-dependent 
Hamiltonian describing the rf-pulses along the x-axis. This sequence implements the 
Hamiltonian of Eq. (1) to the lowest order in an average Hamiltonian theory sense [45]. 

In an inductively detected NMR experiment where we measure quadrature signals 
from the x and y components of the precessing magnetization, the observed signal is 
S(t) = ((a-(t)} = CTr {(T-p(t)}, where a~ = £\<rj = - ia))/2, and C is 

a proportionality constant. The only terms in p(t) that yield a non-zero trace, and 
therefore contribute to S(t), are angular momentum operators such as a^~, which are 
single-spin, so-called "single-quantum coherences" in the language of multiple-quantum 
(MQ) NMR [ 12]. In a standard MQ experiment that is used to characterize the many- 
spin dynamics of a nuclear spin system, the uncorrelated thermal initial states of the 
spins are allowed to evolve under a Hamiltonian such as the DQ Hamiltonian that 
generates the multi-spin dynamics. The multi-spin character of the states is indirectly 
encoded by a collective rotation of the spins and the DQ evolution is then reversed to 
convert the many-spin states back to single spin terms that can be detected. 

Our experiments were performed at room temperature in a 7T vertical bore NMR 
magnet on a single FAp crystal with its c-axis aligned to the external field. The measured 
Ti of the fluorine spins was 300 s (many orders of magnitude larger than the timescales 
explored experimentally here). The length of the n/2 pulse used was w = 1.06 ps 
and A=2.9 /is. The thermal (equilibrium) state of the spins in high field and high 
temperature is highly mixed, and is given by p t h ~ 1 e Si z; where e ~ 10 -5 1. 
The identity is unchanged under unitary transformations, and does not contribute to the 
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Figure 1. Evolution of the thermal initial state under T-L-dq. The blue squares 
correspond to the collective magnetization ■ a\) and the red asterisks to the two- 
spin correlations (J2i( a x <7 y +1 + a y <T x^ 1 ))- The dashed black line is the best fit of the 
observed magnetization to the analytical model described later. The solid red line is 
a guide to the eye. The inset shows the NMR FID of FAp. 

signal S(t). Thus, it is only the deviation of the density operator from the identity that 
gives rise to an observable signal. The constant e then becomes a scaling parameter and 
its value does not affect the details of the experiment (as long as the high temperature 
approximation remains valid). 

Figure 1 shows the observed evolution of the collective magnetization (oc £V a l z - 
blue squares) and the two-spin correlations (oc X^i( a z a y +1 + a y°a; +1 ) _ re d asterisks) under 
"Hdq, starting from p t h- The error bars were estimated from the standard deviation of 
a signal-free region of the time-domain data. The evolution time was incremented by 
increasing the number of cycles from 1 to 40 (30 in the two-spin correlation readout). 

In contrast to typical NMR experiments which involve the DQ Hamiltonian [43], no 
evolution reversal was performed before signal detection in our experiments. In order 
to minimize receiver dead-time effects, a solid echo was used to read out the (single- 
spin) magnetization terms. The two-spin terms were read out using a 7r/4 pulse to 
generate a dipolar echo (with an appropriate phase cycle). The combination of a 7r/4 
pulse followed by evolution under the dipolar Hamiltonian refocuses a portion of the 
two-spin correlations back to single-spin coherences that can then be detected [41]. The 
importance of the single-spin and two-spin operators in driving the transport dynamics 
can be seen from the following equations (the short time dynamics are further discussed 
in Appendix B): 

jfl = y (?\ol + a\a x 2 ) , j t (afa y 2 + a\a x 2 ) = -da z 2 + 3 spin terms. (3) 
As a result, the single-spin and two-spin terms are observed to evolve 90 degrees out 
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Figure 2. Evolution of the end-polarized initial state under Huq- The blue squares 
correspond to J2i a z> while the red asterisks correspond to Yli( a x a y +1 +< T y <J l +1 )- Inset 
shows the 4-pulse experiment on the thermal initial state. The dashed black lines are 
the best fits of the magnetization dynamics to the analytical model. 



of phase with each other during DQ evolution. The inset shows the observed NMR 
Free Induction Decay (FID), which corresponds to the evolution of ^ a l x under "Hoip- 
While the FID decays in about 350 /is, due to the creation of multi-spin correlations 
[18], the magnetization oscillations persist for up to ~ 1.5 ms under Hdq, indicating 
that high-order spin correlations develop quite slowly [46]. 

It is also possible to experimentally prepare initial states of the form puv ~ 
1 — t[o\ + o~ z N ) [47], in which the polarization is localized at the ends of the chain. 
We call this the end-polarized state. Figure 2 shows the DQ evolution observed for pi^. 
Here, we implemented Hdq using the first 4 pulses of the above sequence (C ■ C) to 
achieve better temporal sampling of the signal, though this introduces first-order errors 
in the resulting average Hamiltonian. The alternating sign of the local magnetization 
during spin transport of the end-polarized state results in a rapid attenuation of the 
amplitudes of the observed signal (the 8-pulse version of the experiment is included in 
Appendix A). For comparison, the inset shows the 4-pulse version of the thermal state 
experiment. 

3. Theoretical Analysis and Interpretation 

3.1. Fermionic Model 

Both the NN XX [7] and NN DQ [47, 48] Hamiltonians are well known to be analytically 
solvable in ID, by means of a Jordan- Wigner mapping onto a system of free fermions. 
We invoke this mapping to interpret the experimental results. Starting from an end- 
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polarized initial state pj = 1 — crj, we once again take note of the fact that the identity 
does not evolve, and focus our attention on the deviation term crj. The deviation density 
operator at time t under the DQ Hamiltonian is given by [47, 49]: 

N 



^ = jVTT ^ sm(kj)sm(hj) x (4) 

k,h=l 

a\a h + a[a k ) cos ip kh (t) - i \ a\a\ - a h a k ) smi/; kh (t) 



where a k = J Ylh=i ^' m (kh)ch, Ch = — Yl^i^i) 17 - are canonical fermionic 
operators, and ip k h(t) = 2dt[cos(k) + cos(/i)], and we have assumed open boundary 
conditions on the chain. In order to characterize the evolution of the individual spin 
operators, we can express the evolved state in terms of the Ch operators. This 

yields: 



p J DQ (t) = (-ir 1 { J2 {clc q + c\c p )A hq (t)A hp (t) 

p—q S even 



+ E {P ~ q (°l c l - A^A^t)}, (5) 

p—q 6 odd 

where the time-dependent amplitudes Aj >(? read 

oo 

AM = ^JlmN^dt) ~ i*J 2m * +E (2dt)] 

m=0 

oo 

+ £ i 2m * \i- 5 J 2m ^dt) - i^J 2m ^{2dt)] (6) 

m=l 

with N = N + 1, 5 = q — j, £ = g + j, and J n ( 4 ) being the n-th order Bessel function 
of the first kind. Similarly, for a thermal initial state p t h — 1 — Y2j °z> we ^ n( ^ the 
deviation density operator at time t 

p t D h Q w = - E (4 c o + c l°p) A ^ 2t ) 

p—q S even 

+ E {clc\-c q c p )A M (2t). (7) 

p—q S odd 

Taking the even/odd constraint into account, the density operators in Eqs. (5) and 
(7) are seen to belong to a N(N + 1)/2- dimensional operator subspace which defines the 
Liouville space within which the transport occurs. Mapping back to spins, and assuming 
p > q, we get: 

C p C q = a q a q+l ' ' ' a p-l G p1 C q C P = °q a q+l ' ' ' °p-\°p i I o\ 

c\c P = or+a z q+1 ■ ■ ■ a z p _ x o7 p , c\c q = § (l - cr|) . 
The above quadratic scaling is the same that Fel'dman and coworkers theoretically 
established for the XX Hamiltonian [50] J. While pth is a constant of the motion under 

X It is also interesting to observe that the operators in Eq. (8) are the same as the string operators in 
the spin-spin correlation functions Sxx(<7jP) defined in T. S. Cubitt and J. I. Cirac, Phys. Rev. Lett. 
100, 180406 (2008). 



ID Magnetization transport 



8 



the XX Hamiltonian, evolution of the end-polarized state can be readily expressed as 

pf x (t) = - E iP+q ( c K + c W) -wo. wo- (9) 

3.2. Magnetization Dynamics and Transport Velocity 

We can re-examine the experimental data in Figs. 1 and 2 using the analytical results 
outlined above. The magnitude of the collective spin magnetization is given by 
Sit) = 52p=i ^i,pC0 f° r an end-polarized initial state, and by S(t) = J2p=i A p , p {2t) 
for a thermal initial state. The dashed black lines in the figures are the best fits to these 
expressions, where we have assumed that N is sufficiently large that no boundary effects 
are observed. Thus we only used the m = term in Eq. (6) to calculate Ai tP and A IhP to 
calculate S(t) as shown above. For the end-polarized state we have also assumed that 
on the timescale of the experiment the magnetization propagating from the two ends 
have not had a chance to overlap, and we can thus describe them as two independent 
chains. This allows us to ignore the A^^ term. 

Three fitting parameters were used: a scalar multiplier, the frequency argument 
of the Bessel function, and an additive baseline constant. The baseline constant was 
subtracted from the data shown in the two figures, so that the oscillations are observed 
around zero. For the thermal initial state, the observed signal is seen to damp out at a 
faster rate than expected from the model. This is likely due to the presence of longer- 
range couplings that have been ignored in the NN model, as these long-range couplings 
lead to a leakage out of the restricted Liouville space. The obtained fitting frequencies 
yield d = 8.32 x 10 3 rad/s (8-pulse, thermal state), d = 8.52 x 10 3 rad/s (4-pulse, 
thermal state), d = 8.71 x 10 3 rad/s (8-pulse, end polarized state - see Appendix 
A), and d = 9.56 x 10 3 rad/s (4-pulse, end polarized state), yielding an estimate 
d = 8.78 ± 0.55 x 10 3 rad/s. The estimate is biased by the value obtained in the 
4-pulse experiment on the end-polarized state, where the fitting frequency is seen to 
be too high at longer times (Fig. 2). However, the estimate is still in good agreement 
with the values d = 8.3 x 10 3 rad/s, obtained from observing MQ coherences [ 10], and 
d = 8.17 x 10 3 rad/s, obtained from the known structure of FAp [38]. 

In the thermodynamic limit, the magnitude of the magnetization transported from 
site j — 1 at time t — to site j — n at time t is given by Py° n {t) = Af n (t) 2 , where 

nJ n {2dt) . n ^ 1 J n . 1 {2dt) + J n +i(2dt) 
Al ' n{t) ~ % 2dt ~ % 2 ' { } 

Using properties of the Bessel functions [5 1 ] , it is possible to show that d t Af n = 
— d(A^ n+1 — A™ n _]). If we define a continuous spatial variable z = aj, with a equal 
to the distance between two spins and j being the spin number in the chain, we can 
replace the finite difference with a spatial derivative d t Af n = —2add z A^° n . Taking the 
second derivative with respect to time we thus obtain a wave equation for the transport 
amplitude: 

d 2 t A™ n = [2adfdlA™ n . (11) 
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P^° n (t) also follows the same wave equation with velocity lad. We can also calculate 
the group velocity of the spin system directly from the dispersion relation for the DQ 
Hamiltonian, u(k) = 2d\ cos(fca)|, which results in v g = doj/dk\k=n/2a = lad. Using 
d = 8.78 ± 0.55 x 10 3 rad/s, and a = 3.442 A, we obtain v g « 6.04 ± 0.38 /rni/s. 
This corresponds to a displacement of ~ 9.07 ± 0.57 nm in 1.5 ms or transport across 
N = 26 ± 2 spins. The time taken to travel a distance of n lattice sites is t — n/2d. At 
this time A™ n = i n - l J n {n), and P£ n = {-l) n ~ l J 2 n {n) . For large n, J n {n) ~ ^ (^) 1/3 
[51] and the magnitude of the polarization transported from spin 1 to n scales as 

Pr n (n) ~ n-V\ (12) 



3. 3. Evolution of Multi-spin Correlations 



Figure 3 shows the dynamics of the 1,2,3, and 4-spin correlations for aiV = 9 spin chain 
initialized in the state p\ = 1 — o~\. These curves were generated using Eqs. (3) and (4). 
We can express Ai jQ as 



E 



2mN+q-l 



2mN + q 
dt 



■2mN—q+l 



2mN — q 
dt 



Jo 



2mN-q\^dt) 

z) [ ]. The 



2mN+q 

m=0 \ J m=l 

where we have used the Bessel function identity J v -i(z) + J u +i(z) = ^fJ v 
sequential growth outlined above can easily be visualized in these simulations. The 
figures show that the presence of the boundary (at 2dt m w 5, where t m ~ N/2d 
identifies the so-called "mirror time", see also next section) has a large effect on the 
spin dynamics. In Appendix B we show that additional insight into the dynamics 
is obtained by calculating the short-time evolution, and examining the growth of the 
multi-spin operators. 

It is possible to experimentally characterize the growth of these multi-spin 
correlations. In a MQ experiment [13], we record the coherence orders by performing 
a collective rotation of the state about the z-axis, and observing the resulting phase 
shifts. The first row of Eq. (8) shows the double-quantum coherences created during 
NN DQ evolution, while the second row shows the zero-quantum coherences and the 
polarization states. 

While we cannot directly observe the high-spin correlations in Eq. (8), we can 
use so-called x-basis encoding techniques to characterize the distribution of p — q [44] , 
thereby indirectly probing the growth of these terms. For example, the first term in Eq. 
(8) in the x-basis reads: 



„IE I I 'X- 



+ 0-1 



where 



o-y ± ia\ 



(2i)P-«- 2 

Here, a collective rotation of the system about the x-axis 
results in overlapping binomial distributions of phase factors whose highest order is p—q. 
Higher-order coherences in the x-basis are thus a signature of the presence of multi- 
spin correlations. Figure 4 shows the relative x-basis coherence intensities measured 
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Figure 3. This figure illustrates the evolution of the 1,2,3, and 4-spin correlations 
for a N = 9 spin system initialized in the state p\ = 1 — a\. The amplitudes of 
different product operators are obtained from Eqs. (3) and (4) in the main paper, (a) 
Single spin polarization terms a^. The amplitudes are given by A\ (t). (b) Two 

(2) 

spin correlation terms h n n+1 . These are double quantum terms whose amplitudes 
are given by A 1: „(t)A l ll+1 (i). (c) Three spin correlation terms h^ n+1 n+2 . These are 
zero quantum terms whose amplitudes are given by Ai n (t)Ai n +2(t). (d) Four spin 
correlation terms h!^ n+1 n+2 n +3,- These are double quantum terms whose amplitudes 
are given by Ai^ n {t)Ai tn+ ^{t). Higher order correlations are not shown. 



as a function of the DQ evolution time, starting from the end-polarized state. The 
experiments were performed using a 16-pulse implementation of the DQ Hamiltonian 
[44], and both the number of cycles and the delay A were varied. It can be seen that 
following an initial rapid creation of 3-spin correlations (and concomitant reduction in 
the single spin term), the coherence orders change quite slowly. 



3.4- Mirror Times 

Equations (3)-(7) also help us understand the source of the signal enhancements 
observed near the boundaries of a finite spin chain. The sum over m in Eq. (6) resembles 
a sum over an infinite number of copies of the iV-spin chain (consistent with the periodic 
boundary conditions imposed by the sine transform [49]). At short times, only the 
lower-order Bessel functions contribute, and each "replica" of the chain is independent 
of the others. At longer times, the adjacent copies begin to interfere with each other as 
illustrated in Fig. 5. It is the interference between these terms that is responsible for 
the mirror times observed in our previous work [40], consistent with the interpretation 
in terms of bouncing spin-wave packets and "erratic" dynamics put forward for the XX 
model [50]. In particular, the leading mirror term yields a factor of 4 increase in the 
magnitude of the transferred magnetization. A similar transport behavior has been 
shown for pure states [53]. 
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Figure 4. Growth of x-basis coherences during DQ evolution, starting from an end- 
polarized initial state. Only odd-order coherences are observed due to selection rules. 
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Figure 5. Five copies of a 20-spin chain illustrate the origin of the interference patters 
observed following reflection off the boundary during evolution under the NN XX 
Hamiltonian. The time axis is in units of 1/d. 



Fig. 6 shows that the mirror enhancements obtained under the DQ Hamiltonian 
are also due to interferences generated as a polarization wave-packet reflects off the 
boundary in a finite spin chain. In the case of the DQ Hamiltonian the polarization 
is inverted at every adjacent spin site as it travels down the chain, and the nature of 
interferences at the boundary depends on whether the chain contains an even or odd 
number of spins. 

The above mirror dependence on chain length is also manifested in the multiple 
quantum coherence dynamics of the spin chain. Mirror times were also observed in the 
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Figure 6. Five copies of a spin chain illustrating the boundary interferences for (a) 
N = 20; (b) N = 21 during evolution under the NN DQ Hamiltonian. The time axis 
is in units of 1/d. 

multiple quantum dynamics for thermally polarized spin chains [40]. We can explain 
this as follows. Although significant overlap for waves originating at symmetric spin 
positions j and N + 1 — j can occur at different times depending on j (for example, 
at times tj oc N ± j) there is also large overlap at the time t m oc N, when one of the 
waves has bounced off a boundary and encounters the second wave originating from the 
symmetric spin (see Fig. 7). Since at each time tj there occurs only the overlap for one 
pair, weighted by a factor N, only the overlaps at the mirror time t m oc N are noticeable 
for the thermal initial state [40]. 

4. Summary 

In summary, we have measured the single-spin and two-spin Liouville-space operators 
that drive the (nearly) ballistic transport of magnetization under the DQ Hamiltonian 
at room temperature. We also developed a real-space description of the spin dynamics 
that is able to accurately describe the observed results, by allowing us to characterize 
the growth and coherent dynamics of multi-spin correlations, and to explain the origin 
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Figure 7. Zero quantum coherence intensities Jji(t) + J%+i-ji(t) as a function of 
time in a chain of iV = 31 spins. Left top: j = 1; Left bottom: j = 5; Right top: 
j = 10; Right bottom: j = 15. The black lines are drawn as a guide for the eye. 

of the previously described "mirror times" in finite spin chains. Within the validity of 
a NN approximation, the dynamics of the iV-spin system are seen to be restricted to 
a Liouville space whose size grows only quadratically with N. Based on our model, 
we estimate that the magnetization is coherently transported down the chain with a 
group velocity of 6.04 ±0.38 /im/s. We expect this velocity to be directly independently 
measurable using reciprocal space NMR methods [54] . 
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Appendix A. Evolution of the end-polarized state under "Hog 

Figure Al shows the experimental results obtained using the 8-pulse implementation 
of the DQ Hamiltonian (in contrast to the 4-pulse version shown in Figure 2), when 
the spins are prepared in the end-polarized state piN ~ 1 — e(of + a z N ) in which the 
polarization is localized at the ends of the chain. Fitting the magnetization dynamics 
to the analytical model results in a value of d = 8.7 x 10 3 rad/s. 



ID Magnetization transport 



14 



3000 



2500 




Figure Al. Evolution of the end-polarized initial state under Hdq, showing 
the collective magnetization a\ — blue squares) and the two-spin correlations 
(J2i( a x <J y +1 + < 7 y a x +1 ) — re d asterisks). The dashed black line is the best fit of the 
observed magnetization to the analytical model described later (yielding d = 8.7 x 10 3 
rad/s). The solid red line connecting the red asterisks is a guide to the eye. 



Appendix B. Short time evolution 



Let the initial state of the spin system be p 3 ; = 1 — o{. The state of the spin system at 
a time t later, following evolution under the DQ Hamiltonian is 

t 2 it 3 

Pt = Po-it[7ij3Q,po]-—[H D Q, [Ht>q,Po\]--^- [H dq , [H D q, [H D q,Po]]]±- • • (B.l) 

2 o 

Evaluating the above commutators, we get 

[Hbq, p ] = 2d (hfl lfj + hf} +1 ) , (B.2) 

where hfj = or^a^ — v^aj = | (cfcr| + (rfcrf). A collective rotation of hfj about the 
z-axis by an angle results in phase shifts by angles ±20, indicating a double quantum 
coherence. Similarly, 

[H DQ , [H D q, po]] = - 2d 2 {u]_ x + 2a] + a z j+1 ) 

- 2d 2 (hf2 2 j-ij + 2hf\ j>j+1 + hf j+l j+ ^j , (B.3) 

where hf^ k = afa']a k + o]a\ and 

[H DQ , [H dq , [Hdq, Po]]] = 6d 3 (hf_ %J _ x + 3hf! hj + 3hf] +1 + h$ 1J+2 ) - (B.4) 

where h^./ = °~t°) ~k t + a 7 <T j cr k '7 ^ an< ^ so on - Collective rotations of hfj k and 
^iifc.z about the z-axis by an angle results in phase shifts by angles and ±20 
respectively, indicating these are zero and double quantum coherences. We observe 
the same alternating zero and double-quantum signatures in the commutator expansion 
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as obtained in the fermionic solution (Eq. (3) in the main text). Moreover, here we 
see that the growth of the spin system is necessarily slow. At the level of the second 
commutator, we can infer that the system is 3 times more likely to evolve from a 3-spin 
state to a 2-spin state, than it is to evolve to a 4-spin state. This is one indication of 
why locally polarized states are transported almost ballistically along the spin chain. 
This is schematically illustrated in Fig. Bl. 



1 Z 

H 

2 XY+YX 
It 

3 XZX+YZY 
XZZY+YZZX 



Figure Bl. The successive commutators in Eqs. B1-B4 yield the Liouville space 
operators shown in the figure and illustrate the progressive generation of zero and 
double quantum coherences. 
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